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Abstract 

The memristor is a device whose resistance changes depending on the 
polarity and magnitude of a voltage applied to the device's terminals. 
We design a minimalistic model of a regular network of memristors using 
structurally-dynamic cellular automata. Each cell gets info about states 
of its closest neighbours via incoming links. A link can be one 'conductive' 
or 'non-conductive' states. States of every link are updated depending on 
states of cells the link connects. Every cell of a memristive automaton 
takes three states: resting, excited (analog of positive polarity) and re- 
fractory (analog of negative polarity) . A cell updates its state depending 
on states of its closest neighbours which are connected to the cell via 'con- 
ductive' links. We study behaviour of memristive automata in response 
to point-wise and spatially extended perturbations, structure of localised 
excitations coupled with topological defects, interfacial mobile excitations 
and growth of information pathways. 

Keywords: memristor, cellular automaton, excitable medium 

1 Introduction 

The memristor (a passive resistor with memory) is a device whose resistance 
changes depending on the polarity and magnitude of a voltage applied to the 
device's terminals and the duration of this voltage's application. Its existence 
was theoretically postulated by Leon Chua in 1971 based on symmetry in in- 
tegral variations of Ohms laws [H |5J IE] • The memristor is characterised by a 
non-linear relationship between the charge and the flux; this relationship can 
be generalised to any two-terminal device in which resistance depends on the 



internal state of the system [5]. The memristor cannot be implemented us- 
ing the three other passive circuit elements — resistor, capacitor and inductor 
therefore the memristor is an atomic element of electronic circuitry [4| [Sj, E] ■ 
Using memristors one can achieve circuit functionalities that it is not possible 
to establish with resistors, capacitors and inductors, therefore the memristor is 
of great pragmatic usefulness. The first experimental prototypes of memristors 
are reported in [191 [7J 120) . Potential unique applications of memristors are in 
spintronic devices, ultra-dense information storage, neuromorphic circuits, and 
programmable electronics |17j . 

Despite explosive growth of results in memristor studies there is still a few 
(if any) findings on phenomenology of spatially extended non-linear media with 
hundreds of thousands of locally connected memristors. We attempt to fill 
the gap and develop a minimalistic model of a discrete memristive medium. 
Structurally-dynamic (also called topological) cellular automata [11], [8] seem 
to be an ideal substrate to imitate discrete memristive medium. A cellular 
automaton is structurally-dynamic when links between cells can be removed 
and reinstated depending on states of cells these links connect. Strcturally- 
dynamic automata are now proven tools to simulate physical and chemical dis- 
crete spaces [TUGS 0331 USE] an d graph-rewriting media [18]; see overview in [j"2] . 

We must highlight that simulation of cellular automata in networks of mem- 
ristors is discussed in full details in [TJ] . Itoh-Chua memristor cellular automata 
are automata made of memristors. Memristive cellular automata studied in 
present paper are cellular automata which exhibit, or rather roughly imitate, 
certain memristive properties but otherwise are classical excitable structurally- 
dynamic cellular automata. 

The paper is structured as follows. Memristive automata are defined in 
Sect. [2] In Sect. [3] we analyse space-time dynamics of the automata in response 
to external perturbations. Anatomy of stationary oscillating localizations is 
presented in Sect. [4] We analyse structure of disordered excitation worms prop- 
agating at the boundary between disorganised and ordered link domains in 



Sect. 
Sect. 



Collision-based approach to layout of conductive pathways is given in 
Some ideas of further studies are outlined in Sect. [7J 



2 Memristive automaton 

Definition 1 A memristive automaton is a structurally- dynamic excitable cel- 
lular automaton where a link connecting two cells is removed or added if one of 
the cells is in excited state and another cell is in refractory state. 

A cellular automaton A is an orthogonal array of uniform finite-state ma- 
chines, or cells. Each cell takes finite number of states and updates its states 
in discrete time depending on states of its closest neighbours. All cells update 
their states simultaneously by the same rule. We consider eight-cell neigh- 
bourhood and three cell-states: resting o, excited +, and refractory — . Let 
u(x) — {y : \x — j/|loo = 1} be a neighbourhood of cell x. A cell x has a set of 
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incoming links {l xy : y £ u(x)} which take states and 1. A link l xy is a link 
of excitation transfer from cell y to cell x. A link in state is considered to be 
high- resistant, or non- conductive, and link in state 1 low resistant, or conduc- 
tive. A link-state l xy is updated depending on states of cells x and y at time 
step t: l xy = f{x t ,y t ). Resting state gives little indication of cell's previous 
history, therefore we will consider not resting cells contributing to a link state 
updates. When cells x and y are in the same state (bother cells are in state + or 
both are in state — ) no 'current' can flow between the cells, therefore scenarios 
x 1 = y l are not taken into account. Thus we assume that the only situations 
when x , y G {+, — } and x ^ y may lead to changes in links conductivity: 

l x + y x — { a , x l = + and if = -6, x* = - and y* = +l xy , otherwise (1) 

where a^b and a,b <E {0,1}. Thus we consider two types of automata A ab : 
A 01 and A 10 . 

A resting cell excites (x* = o — > x = + transition) depending on number of 
excited neighbours. 

There are two ways to calculate a weighted sum of number of excited neigh- 
bours: 

!• s + = J2yeu( x )X(y t ,+) 

where +) = 1 if y* = + and x(y*; +) = otherwise. 

Thus, we consider two types of memristive automata. In automaton A\ a 
resting cell excites if cr + > 0. In automaton A2 a resting cell excites if £+ > 1 
or cr + > 0. 

Note 1 Automaton A\ is a pre-memristive cellular automaton because it imi- 
tates only polarity (links update |Ip ) not voltage; automaton Ai is a memristive 
cellular automaton because it imitates both polarity (links update ) and cur- 
rent intensity (use o/E + and <J+). 

A polarity of a current is imitated by excitable cellular automaton using 
exctited + and — refractory states. If a cell x is in state — and a cell y is in 
state + then cell y symbolises an anode, and cell x a cathode. And we can say 
that a current flows from y to x. Indeed, such an abstraction is at the edge of 
physical reality, however this is the only way to develop a minimal discrete model 
of a memristive network. In automaton Ai the condition E + > 1 symbolises 
propagation of a high intensity current along all links, including links non- 
conductive for a low intensity current. This high intensity current in Ai resets 
conductivity of the links and also states of cells. The condition cr+ > reflects 
propagation of a low intensity current along conductive links. The current of 
low intensity does not affect states of links but only states of cells. 
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2.1 Experiments 



We experiment with 300 x 300 cell automaton arrays, with non-periodic ab- 
sorbing boundaries. We conduct experiments for two initial conditions on links' 
'conductivity' — 

• Li-condition: all links are conductive (for every cell x and its neighbour 
y l° y = 1), and 

• Lo-condition: all links are non-conductive (for every cell x and its neigh- 
bour y 1% = 0). 

While testing automata's response to external excitation we use point-wise 
and spatially extended stimulations. By point-wise stimulation we mean exci- 
tation of a single cell (A\) or a couple of cells (A2) of resting automata. These 
are minimal excitations to start propagating activity. 

Let D-disc be a set of cells which lie at distance not more than r from the ar- 
ray centre (n/2,n/2), n = 300, r = ri/3). When undertaking D-stimulation we 
assign excited state + to a cell of D with probability 0.05. In some cases we ap- 
ply D-stimulation twice as follows. An automaton starts in L\- or ^-condition, 
we apply /^-stimulation first time (we call it ^-excitation) and wait till exci- 
tation waves propagate beyond boundaries of the array or a quasi-stationary 
structure is formed. After this transient period we apply D-stimulation again 
(^2-excitation) without resetting links states. 

Space-time dynamics of automata is illustrated by configurations of excita- 
tions and dynamics of link conductivity is shown either explicitly by arrows (in 
small configuration) or via grey-scale representation of cells' in-degrees: cell x 
is represented by pixel with grey- value 32 • X)j/eu(xV Despite not representing 
exact configuration of local links, in-degrees give us a rough indicator of spatial 
distribution of conductivity in the medium. The higher is the in-degree at a 
given point, the higher is the conductivity at this point. 

3 Phenomenology 

Finding 1 A point-wise stimulation of automaton Ai leads to a persistent ex- 
citation, while automaton A\ returns to a resting state. 

A single-cell excitation of resting automaton A\ (Fig. [I]) or two-cell excitation 
of resting automaton A% (Fig. [2]) in L\ initial conditions lead to formation of a 
'classical' excitation wave-front. The wave-front propagates omni-directionally 
away from the initial perturbation site and updates states of links it is passing 
through. 

Links leading from cells to the neighbours they excited are made non-conductive 
in development of A® 1 (Figs. [l] and [2]); or we can say that links corresponding to 
normal vectors of propagating wave-front are made non-conductive. Cell excited 
at time step t = 1 becomes isolated. The situation is similar in development 
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(a) t = 1 (b) t = 2 (c) t = 3 



\XT>dXTxTxXXX Xr>dVlXX^T><XXX XXTxTxTxTxtxl/ 




/IXIXIXJAIXLXJX XlxyuXlXlXlXX / I^IXI^4A1XJXJ \ 

(d) t = 4 (e) t = 5 (f) t = 6 



Figure 1: Snapshots of excitation and links dynamics in automaton A® . In 
initial configuration (t = 1) one cell is excited others are resting. Arrow from 
cell x to cell y means l xy = 1. Excited cells are shown by red discs with sign 
'+', refractory cells by blue discs with sign resting cells are blank. Links at 
the edges of cellular array corresponds to an absorbing boundary cells, which 
are always in resting state independently on states of their neighbours. 



of automata A\° and A^° with the only difference that links connecting cells 
which are excited at time step t to cells they have been excited by are removed. 

In summary, in automata A® 1 links associated with forward propagation of 
perturbation are made non-conductive, and in A} links associated with back- 
ward propagation are made non-conductive. Excitation wave-front travelling 
from a single stimulation site forms a domain of co-aligned links. Excitation 
waves initiated in different cells collide and merge. Boundaries between do- 
mains formed by different fronts are represented by distinctive configurations of 
links (Fig. [3]). 

Automaton A\ is non-excitable in Lo-conditions because no excitation can 
propagate along non-conductive links. An outcome of two-site excitation of 
resting automaton Ai in Lo-condition depends on configuration of the initial 
excitation. 

Let two diagonal neighbours (north-west and south-east) be excited (Fig. [4^,) 
at time step t. These two cells transfer excitation to their two neighbours (north- 
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/!AIAIAIAIAP<JAI\ / !XlAtAlAlAP<4A] \ 
(a) t = 1 (b) £ = 2 




/!AIAIAIAP<4AIAI\ / !Xl^JAlAP<a><JX] 
(c) t = 3 (d) t = 4 



Figure 2: Snapshots of excitation and links dynamics in automaton A^ 1 - In 
initial configuration two cells (minimum size of non-dying excitation for A2) are 
excited others are resting. Arrow from cell x to cell y means l xy = 1. Excited 
cells are shown by red discs with sign '+', refractory cells by blue discs with 
sign resting cells are blank 

east and south-west) as shown in Fig. |4Jd. Only links from north-west cell to 
north-east and south-west, and south-east to south-west and north-east are 
formed (Fig. [3J:) and excitation becomes extinguished (Fig. |2Ji). 

Let excited cells be neighbours at the same row of cells (Fig.p)k). If they are 
excited at time step t then their north and south neighbours are excited due to 
condition S + > 1 taking place (Fig.^). The localised (two-cell size) excitations 
propagate north and south (Fig.[5j;) and make links pointing backwards (towards 
source of initial excitation) conductive. At the second iteration cell lying east 
and west of initially perturbed cells becomes excited (Fig.[5p). By that initially 
perturbed cells return to resting state (Fig. [5ji) and thus they become excited 
again. A growing pattern of recurrent excitation fills the lattice (Fig. [5jic). 

Finding 2 Repeated stimulation of memristive automata in a spatially- extended 
domain leads to formation of either disorganised activity domain emitting target 
waves of excitation or a sparse configuration of stationary oscillating localiza- 
tions. 
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(a) t = 1 (b) t = 2 (c) t = 3 




(d) t = 4 (e) t = 5 



Figure 3: Snapshots of excitation and links dynamics in automaton A® . In 
initial configuration several cells (those in '+'-state at t = 1) are excited others 
are resting. Arrow from cell x to cell y means l xy = 1. 
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(a) t = 1 (b) t = 2 (c) t = 3 (d) t = 4 

Figure 4: Snapshots of excitation and links dynamics in automaton : Lq- 
start. In initial configuration two cells (minimum size of non-dying excitation 
in A) arc excited others are resting. Arrow from cell x to cell y means l xy = 1. 



Examples of excitation dynamics in response to E\- and ^-excitations are 
shown in Fig. [6] Each singular perturbation in the first stimulation of automata 
being in Lx-condition leads to propagation of excitation. The waves merge into 
a single wave and disappear beyond lattice boundary. Merging of waves outside 



D into a single wave is reflected in cells' in-degree distribution in Fig. 19 
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(d) t = 4 (e) t = 5 

Figure 5: Snapshots of excitation and links dynamics in automaton A2 , Lq- 
start. In initial configuration two cells (minimum size of non-dying excitation) 
are excited others are resting. Arrow from cell x to cell y means l xy = 1. 



Automata' behaviour become different after second D-stimulation (i?2-excitation). 
They all exhibit quasi-chaotic dynamic inside boundaries of disc D (shown 
by circle in Fig. [6|); cell in-degrees' distributions of the automata are similar 



(Fig. 18). However the excitation dynamics in automaton is reduced (after 
long transient period) to stationary oscillating localizations while in automata 
.A} , A® 1 , and A^ excitations outside boundaries of initial stimulation merge 
into target waves (Fig. [6] and Fig. 19). Oscillating localizations developed in 
A® 1 stay inside the boundaries of D. 

Random excitation is extinguished immediately in automata A\ being in 
Lo -con dition of total non-conductivity. When automaton A2 in Lo _con( lition 
is excited, the quasi-random excitation activity persists inside boundaries of D 
while omni-dircctional waves are formed outside D. Patterns of activity are not 
changed significantly after second random perturbation, i^i-excitation (Fig. [6]). 

If there are both excited and refractory states in the external stimulation 
domain the developments are almost the same but A\° and A^ shows persistent 
excitation activity already at the first stimulation. 
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Figure 6: Example of automaton behaviour in response to /^-stimulation. 
Boundaries of disc D are shown as circles. Topology of links is represented 
by grey-levels of cells' in-degrees. 
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4 Oscillating localisations 



Finding 3 E2- excitation of A® 1 leads to formation of excitation wave-fragments 
trapped in a structurally defined domains. 

^-excitation of A® 1 leads to formation of sparsely distributed localised oscil- 
lating excitations, or oscillators (Fig. [7]). An oscillating localisation (oscillator) 
usually consist of one or two mobile localizations which shuffle inside a small 
compact domain of the automaton array. This micro-wave is updating states 
of links and thus influencing its own behaviour, A® 1 , a ft er ^-stimulation and 
formation of oscillating localizations, is characterised by a smooth balanced dis- 



tribution of cells in-degrees (Fig. 18 table entry L±, E2). 

Examples of two most commonly found oscillators — 0\ and O2 — are 
shown in Figs. [§] and [9j and their characteristics in Fig. [lO] Both oscillators 
have exactly the same minimum and maximum masses (measured as a sum of 
cells in excited and refractory states). Oscillator O2 has much longer period 
than 0\ and larger maximum density. Oscillator 0\ spans large space during 
its transformation cycle, it occupies a sub-array of 6 x 6 cells when in its largest 
form. 

In its minimal form 0\ consists of two cells: one cell is in state + another 
in state — (Fig. [8j t — 0). At next two steps of Oi's transformations a small 
localised excitation is formed. It propagates north (Fig.[8j t = 1, 2,3). At fourth 
step of oscillator Oi's transformation the travelling localised excitation splits 
into two excitations: one travels north-west, another south (Fig.|8j t = 4,5). The 
localisation travelling south returns to exact position of the first step excitation, 
just with swapped excited and refractory states (Fig. [8j compare t — 6 with 
t = 0) and then repeats a cycle of transformations t = 1,2 (compare t = 1 
with t = 7 an t = 2 with t = 8). The excitation travelling north-west become 
extinguished (t = 8,9). 

A couple ( — h) is a minimal configuration of oscillator O2 (Fig. |9j t = 0). 
When the localisation starts it is development in configuration ( — h), and config- 
uration of links' conductivity as shown in Fig.[9j it is transformed into excitation 
wave-fragment propagating east and north-east (Fig. [9j t = 1,2,3). At fourth 
step of oscillator's transformation the wave-fragment shrinks and by step t = 6 
the oscillator repeats its original state ( — h) yet rotated by 90° clockwise. New 
excitation wave-fragment emerges and propagates west and south-west (Fig. [9j 
t = 7... 10). The fragment contracts to configuration (H — ) by step t = 12. 
Then development and transformation of excitation wave-fragments is repeated 
(Fig. |9j t = 13, • • • , 23). The localisation returns to its original configuration 
(-+) at t = 26. 

Most oscillating localizations observed in experiments with have very 
long periods. This is because the oscillators' behaviour is determined not only 
by cell-state transitions rules, as in classical cellular automata, but also by 
topology of links modified by repeated random stimulation and dynamics of the 
links affected by oscillators themselves. 
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(a) 




(c) 

Figure 7: Formation of localised excitations in A® 1 , Li-start. A localised exci- 
tation occupies a fixed size domain, the excitation may change its size period- 
ically but never expands more then the certain fixed size, (a) configuration of 
automaton after i? 2 -excitation, boundary of D is shown by black circle; (b) the 
same automaton after transient period, only few oscillating localizations sustain, 
(c) links' conductivity presented via grey-values of cells' in-degrees. 
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t=12 t=13 t=14 t=15 




t=16 t=17 t=18 t=19 




t=20 t=21 t=22 t=23 




t=24 t=25 t=26 



Figure 9: Example of oscillator O2 in .A^ 1 - 
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12 
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36 


9 


Minimum density 
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Figure 10: Characteristics of commonly found oscillators 0\ and O2 in A® 1 . 
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5 Dynamics of excitation on interfaces 



Finding 4 Automaton A® exhibits localizations propagating along the bound- 
ary of disordered and ordered conductivity domains. 

As previously discussed, after I?2-excitation of A® 1 the automaton becomes 
separated into two domains. One domain has a disorganised appearance, it lies 
inside boundaries of D. This domain is characterised by quasi-random discoidal 
configuration of cell in-degrees (Fig. [7};). Second domain, cells lying outside 
-D, represents ordered configurations of links (Fig. [7j;) . The links are co-alined 
during target wave propagation outward D. See Figs. [18] and [19] for distribution 
of cell in-degrees in disordered and ordered domaines. 

In automaton A\ l a link from cell y to cell x becomes non-conductive at 
time step t + 1 if x 1 = + and y* = — . Thus target waves travel ouside of D 
form non-conductive pathways. If we stimulate automaton in these domains of 
ordered links no distal propagation of excitation will be observed, however, a 
worm of disordered excitations is formed at the boundary between disordered 
and ordered configurations of links. 

When configuration of oscillating localizations becomes stationary (in a sense 
that the localizations do not travel outside their fixed domains) we apply a few- 
cells wide excitation at the boundary of D (Fig.[TlJi), north-north- west localised 
excitation, which exceeds in size existing oscillating localizations). An extended 
cluster, or a worm (a group of excited and refractory states with a preferable 
direction of growth), of excited and refractory cells is formed (Fig. |TlJ)) . It 
expands along the boundary between domains. Two heads of the emerging 
worm — one head moves clockwise and another anticlockwise — are clearly 



visible in (Fig. 11). At some point the worm undergoes sub-division (Fig. 11 



and it breaks into two independent worms (Fig. Hi). Each new worm has an 
extended head and gradually thinning tail. Thickness of a worm, measured in 
a number of non-resting cells along the worm body, is shown in Fig. 

The worms move along boundary of D until they collide (Fig. [IT 
nihilate in the result of the collision (Fig. [TT| i) . The configuration of oscillating 
localizations then returns to its orig inal state (Fig. [El}). 
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(a) t (b) t + 200 



f s * 

r - ,.i c 1 - • ..; 



(c) i + 300 (d) t + 500 



r 



(e) t + 600 (f ) t + 800 

Figure 11: Propagation of localised excitations at the interface between disor- 
dered and ordered conductivity domains, A® 1 automaton, Li-start. 



15 





a 



3jt/2 
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Figure 12: Illustrative profile of a worm excitation propagating anti-clockwise 
in Fig. |11| 77(a) is a normalised number of cells in non-resting states lying at 
distance [n/3 — 10, n/3 + 10] and angle a from the centre of disc D. 
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6 Building conductive pathways 



Finding 5 By exploring collisions between excitation wave- fragments travelling 
in A2 one can built information transmission pathways in an initially non- 
conductive medium. Routing primitives realised include signal splitting, signal 
echoing and signal turning. 
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(d) t = 1 (c) t = 1 (f) t = 1 

Figure 13: Snapshots of excitation and links dynamics in automaton , Lo- 
st art. In initial configuration two cells are excited, two cells are in refractory 
state, others are resting. Arrow from cell x to cell y means l xy = 1. 

Localizations travelling in A2, Lo-condition, can form pathways conductive 
for low-strength excitations. For example, a localisation of two excited and two 
refractory states propagates in Ai in the direction of its excited 'head'. The 
localisation forms a chain of conductive links oriented opposite to the localisa- 
tion's velocity vector in case of A%, and in the direction of the localisation's 



propagation in case of A\ (Fig. 13 ) 



A layout of conductive pathways, or 'wires', is determined by outcomes of 
collisions between path-laying particles. A detailed example of such collision- 



determined pathway building is shown in Fig. 14 Two travelling localizations 
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Figure 14: Snapshots of excitation and link dynamics in automaton , L$- 
start. In initial configuration two 2 + -particles are introduced in the medium. 
Four cells are excited, four cells are in refractory state, others are resting. Arrow 
from cell x to cell y means l xy = 1. 



2 + -particlea^[ are initiated in A^ facing each other with their excited heads 
(Fig. [Mfc ,) . One particle propagates south another north (Fig. 14 1) . When the 
particles collide they undergo elastic-like collision, in the result of which two 



1 2+-particles are localizations consisting of two excited and two refractory states, which 
move along rows or columns of an excitable cellular array [2] 
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2 + -particles are formed: one travels east another west (Fig. 14;— f). When a 
weak excitation rule — a cell is excited if at least one neighbour is excited and 
no links are updated — is imposed on the automaton the pathways formed by 
these two colliding particles becomes selective. If an excitation is initiated in 
the southern or northern channel the excitation propagates till cross-junction 
and then branches into eastern and northern channels. 

Pathways built by two 2 + -particles undergoing head-on collision being in 
different phases (odd and even distance between their start positions) and lateral 
offsets are shown in Fig. 15 Most T-bone collisions between 2 + -particles lead to 
formation of omni-directionaly growing excitation patterns. Only in situations 
when one particle hits a tail of another particle no uncontrollable growth occurs. 
The particle hitting a tail of another particle extinguishes and the other particle 
continues its journey undisturbed. 

There are primitives of information routing implcmcntablc in collisions be- 
tween 2 + -particles (Fig. 16 1. T-branching, or signal splitting, (Fig. 16 1) is built 
by particles colliding with zero lateral shift and even number of cells between 
their initial positions (e.g. Fig. 15 entry (even, 0)). When signal travelling 
north reaches cross-junctions it splits into two signals - one travel west and 
another travels east; no signal continues straight propagation across junction. 

Echo primitive (Fig. 16 3) is constructed in a head-on collision between 2 + - 
particles with lateral shift two or three cells (see (e.g. Fig. [l5j entries (odd, 
2), (odd, 3), (even, 2), (even, 3)). The echo primitive consists of two anti- 
aligned (e.g. one propagates information northward and another southward) 
parallel information pathways. There is a bridge between the pathways. When 
signal propagating along one pathway reaches the bridge, it splits, daughter 
signal enters second pathway and propagates in the direction of mother signal's 
origination. In (Fig. |16[ d) mother-signal propagates north and daughter-signal 
south. 

Turn primitive (Fig. 
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is implemented in T-bone collision between two 
2 + -particles. A signal generated at either of the pathways propagates in the 
direction of 2 + -particle which trajectory was undisturbed during collision. 



7 Discussion 

We designed a minimalistic model of a two-dimensional discrete memristive 
medium. Every site of such medium takes triple states, and a binary conduc- 
tivity of links is updated depending on states of sites the links connect. The 
model is a hybrid between classical excitable cellular automata [H] and classical 
structurally-dynamic cellular automata |11) . A memristive automaton with bi- 
nary cell-states would give us even more elegant model however by using binary 
cell-states we could not easily detect source and sink of simulated 'currents'. 
Excitable cellular automata provide us with all necessary tools to imitate cur- 
rent polarity and to control local conductivity. From topology of excitation 
wave-fronts and wave-fragments we can even reconstruct relative location of a 
source of initiated current. 
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We defined two type of memristive cellular automata and characterised their 
space-time dynamics in response to point-wise and spatially extended pertur- 
bations. We classified several regimes of automata excitation activity, and pro- 
vided detailed accounts of most common types of oscillating localizations. We 
did not undertake any systematic search for minimal oscillators though but just 
exemplified two most commonly found after random spatially-extended stimu- 
lation. Exhaustive search for all possible localised oscillations could be a topic 
of further studies. 

With regards to formation of conductive pathways just few possible versions 
amongst many implementable were discussed in the papers. Opportunities to 
grow 'wires' in memristive automata are virtually unlimited. For example, in 
A2 1 , Lo-start, after first D-stimulation (Fig. 17 1 generators of spiral and target 
waves are formed inside D. Boundaries between the generators (they provide a 
partial approximation of a discrete Voronoi diagram over centres of the genera- 
tors) are comprised of cells with high in-degrees. Such chains of high in-degree 
cells can play a role of conductive pathways even if we increase excitation thresh- 
old of the medium. 
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(a) 



(b) 



(c) 



Figure 16: Types of information transfer pathways laid out by interacting trav- 
elling localisations: (a) T-branching, (b) echo, (c) turning. 
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Appendix 



Distributions of cell in-degrees for automata A\ and A2 after E\- and E2- 
excitations are shown in Figs. 18 and [19] Calculations are done separately 



for cells lyings inside stimulation disc D (Fig. 18) and outside the disc (Fig. 19). 
Each chart represents distribution of cell in-degrees, where horizontal axis is a 
number of incoming links and vertical axis is a ratio of cells with given in-degree 
to a total number of cells in the analysed domain. 
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Figure 19: Distributions of cells' in-degrees outside disc D. 
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